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SAFETY  STOCK  ALLOCATIONS 


In  Working  Memorandum  No.  128  a  discussion  is  given  of  the  problem  of 
finding  the  optimal  exponent  T  when  the  safety  stock  is  taken  to  be  proportional 
to  the  T  power  of  the  mean  demand.  A  solution  is  obtained  under  the  assumption 
that  the  distribution  of  the  positive-excess  demand  (the  demand  in  excess  of  the 
mean  demand  in  a  lead  time)  is  a  power  function,  with  exponent  approximately  -  1. 

There  are  two  conditions  on  a  distribution  function  that  are  not  met 
by  this  assumption.  First,  the  function  must  be  no  greater  than  1.0;  in  fact,  in 
this  case  it  must  be  no  greater  than  Fq,  the  probability  that  the  demand  exceeds 
the  mean  in  a  lead  time.  Second,  the  expected  value  of  the  square  must  exist  and 
be  finite,  since  we  want  to  assume  that  the  standard  deviation  of  demand  in  a  lead 
time  is  a  xb  ,  where  x  is  the  mean  yearly  demand.  Of  course,  we  can  modify  the 
power  function  (by  truncation  or  otherwise)  to  meet  these  requirements,  but  the 
problem  then  becomes  more  difficult  and  the  solution  more  complex.  The  solution 
is  so  complex, for  any  admissible  distribution  function  that  I  have  tried  (and 
this  includes  several  modifications  of  the  power  function),  that  it  is  not  possible 
at  this  time  to  say  anything  about  the  form  of  the  solution,  except  Jt^at^ls  an 
Important  parameter,  and  probably  the  most  important  one. 
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The  last  line  on  page  4  should  read: 
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This  Is  obtained  by  using  the  Integral  approximation 
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I  J  Jx  F(k(x))g(x)dx 

o 


for  the  sum,  where  F  is  the  power  function,  k(x)  Is  given  by  the  last  line  on 
page  3,  and  g(x)  is  the  log-normal  density  function  for  the  distribution  of  the 
yearly  demand  among  items. 

The  Power  Function 

Let  u  be  the  demand  in  a  lead  time  and  let 


m  •  m(x) 
a  “  o(x)  -  a  x 

be  respectively  the  mean  demand  in  a  lead  time  and  the  standard  deviation  from 
the  mean  for  an  item  having  yearly  mean  demand  x.  The  reorder  point  for  this 
item  is 


r  ■  m  +  X  x 

Since  \  is  positive  and  x  is  non-negative,  a  stockout  can't  occur  when  u  ^m. 
Hence,  only  the  part  of  the  distribution  in  excess  of  m  is  of  interest. 

To  avoid  writing  F  ,  the  probability  that  the  demand  exceeds  the  mean, 
explicitly  as  a  factor,  consider  the  conditional  upper  distribution  function  H(y)  , 
given  that  the  demand  exceeds  the  mean,  where 

y  "  u  -  m. 
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Then,  for  y  ;>  0,  H(y)  is  monotonic  non- increasing  with 

H(0)  -  1,  H(  »  )  -  0 


and 

2  P  00  2 

(1)  a  m  y  dH(y) 

^  o 

Here,  a  is  the  standard  deviation  from  y  “  0,  since  this  value  of  y  corresponds 
to  u  ■  m;  and  it  is  taken  only  on  the  positive  side  of  the  mean,  since  only  these 
values  of  u  are  pertinent. 

Assume  that  H(y)  is  a  function  F(k)  of 

k  -  y/o 


By  this  is  meant  that  when  y  is  replaced  by  ka  the  distribution  function  depends 
only  on  k,  not  on  o.  Then  F(G)  “  1,  F(  »  )  ■  0,  and  (1)  becomes 


(2) 


dF  (k) 


Otherwise,  F  is  a  function  of  both  k  and  a. 

In  Working  Memorandum  No.  128  it  is  assumed  that 

F(k)  -  \k£  ,  e  <0, 
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where  \  and  c  are  independent  of  a.  This  is  not  an  admissible  function,  since 
it  is  greater  than  1.0  for  k  <  k^  ,  where  k^e  ■  1. 

By  truncation  we  get  the  modified  function 


(3) 

which  has 

(4) 


F(k)  - 


the  density  function 


r°  • 

f(k)  »  | 

-\€k 


k  <kl 
k  ^  k1 


k  <  k 

€  1  ,  k>  k 


1 


The  integral  in  (2)  exists  only  if  €  <-2,  unless  we  want  to  truncate  also  at 
some  value  of  k  ■  k2>  The  function  in  (3)  with  £  <-2  is  an  admissible  distribution 
function.  However,  it  is  not  acceptable  as  a  reasonable  distribution,  since  it  is 
not  likely  that  there  will  be  no  demands  near  the  mean,  as  Indicated  by  (4).  In 
fact,  one  would  expect  that  this  region  would  have  highest  density.  Also,  the 
integral  approximation  for  now  would  be  much  more  complicated,  since  it  would 
be  a  sum  of  two  Integrals,  each  of  which  Involves  an  error  Integral. 

A  Modified  Power  Function 

An  acceptable  modification  of  the  power  function  is 


f(k)  = 


f£i  ’  0^k*k2 

-  ek6  £  ,  k  >  l<2 


where 


c  k 


€-1 


e  <-2 . 
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Then 


F(k) 


1  -  fLk 
^  k€ 


0  <  k  <k0 
k  >  k2 


and 


1  -  fjk2  -  ^ 


-  f]k2/e 


or 


f^d-l/e  )  -  1. 


Another  condition  on  the  parameters  is  given  by  (2). 
Now 


k  ■  k(x)  ■  c  x^  k 


where 


c  »  c(r  )  -  (Is/aN)(x)  ‘T  n"  r(  r_1) 
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In  terms  of  x  the  distribution  function  is 


F(k(x)) 


1  -  fjC  X 


r-b 


Xce  x«<r-b) 


0  <c  x  <  x^ 
X>  X2 


where 


k2  -  c(r)x2r"b 


The  integral  approximation  of  becomes 


S  *  * 

1  K 


x1^2  (1  -  f^cx^  b)  g(x)  dx  + 


Y  e  r  1/2+  e  (r-b)  .  x  . 

\c  x  '  g(x)dx 


If  g(x)  is  the  log-normal  density  function,  can  be  expressed  in  closed  form 
in  terms  of  the  error  integral.  Similarly,  the  expected  amount  S2  of  the  shortage 
can  be  expressed  in  closed  form.  But  the  forms  will  be  so  complicated  that  there 
is  little  hope  of  obtaining  the  minimums  analytically.  A  solution  for  any  given 
set  of  parameters  could  be  found  by  approximation. 


Rectangular  Distribution 

The  only  admissible  distribution  F(k)  that  I  have  been  able  to  treat  in 
a  semi-analytical  manner,  when  g(x)  is  the  log-normal  density  function,  is  the 


Arthur  JD.liitU.Jnt. 


WORKING  MEMORANDUM  NO.  136  CASE:  C-63687  DATE:  June  6,  1963  PAGE: 


7 


rectangular  distribution, 

1  -  k/  >JT 

F(k)  -  | 

^  0 


o  <k  ^  JT 
k  >jr 


The  factor  is  needed  to  satisfy  (2);  that  is,  starting  with  f(k)  ■  1/k^  for 
0  k  ^k^,  and  0  elsewhere,  equation  (2)  yields  k^  ■  3 .  Then 


S 


N 

K 


r-b 


C  X 


g(x)  dx  , 


where  (x^.x^)  is  the  interval  of  x  for  which 


o  <cxr'b  <?jr. 


If  T>  b,  x^  ■  0  and  c  x^*  -  JJ".  If  T<b,  c  x^”**  -  JT  ,  and  x2  -  *  . 

For  T>b  the  expression  for  becomes 


Sj  -  (N/K)ic  1/2  n'1/4 


$<  7  log/))  -  n(r‘b)  /°'<r"b>$  (  i  iog/»-8(r-b)) 


where 


s  ■ 


J  2  log 
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/°w  X2/x  -  (l/x)(  JT/c)1^”1^  ,  c  -  c(r) 


\  (u)  -  (1/fU) 


dt 


For  T<b,  the  (p  functions  are  replaced  by  their  complements. 
The  derivative  is 


dsl  . 

r 

1+c ^  J  2 «  exp(z2/2)  J(z)]  ,  1T>b 

dr 

^  A 

j-l+Cj^J  2rc  exp(z2/2)(l- J  (z))  J  ,  Tcb, 

where  A  >  0  and 

(5) 

z  ■  c 

x+c2/w  -  w/2, 

(6) 

w  -  s(r-b)  , 

cL  -  s(b-l/2) 

c2  ■  log(Na  vfVlg)  +  b  log  x+  b(b-l)  log  ji 


We  can  find  the  value  of  z  that  minimizes  and  then  substitute  in  (5)  and  (6) 
to  find  T. 

The  results  arr  as  follows: 

(a)  If  c,<0,  there  is  no  finite  solution.  In  this  case,  I  is  so  large  that 

Z  5 

it  does  not  impose  an  essential  restriction. 


Arthur  2D.Ulttle.3nr. 


WORKING  MEMORANDUM  NO.  136  CASE:  C-63687  DATE:  June  6,  1963 


PAGE:  9 


(b)  If  c2  -  0,  T  - 


(c)  If  c2>0,  the  form  of  the  solution  depends  on  the  sign  of 


(1)  If  ■  0,  that  is,  b  ■  1/2, 
then  T  ■  b  ■  1/2 


(2)  If  C],  J*  0 


b-  7  (  Jj2  +  2  c2  -  J)  ,  b  >  1/2 


{  ; 


b+  £  (  J  J2  +  2  c2  -  J)  ,  b  <  1/2 


where 


J  -  J(  |  cl  I  )  -  z  -  cx 


is  graphed  in  Figure  1.  To  find  T,  compute  the  value  of  |  |  *  |b  -  1/2 

read  the  value  of  J  from  Figure  1,  and  substitute  in  the  above  formulas. 
If  4  c2  <  ,  f  i s  approximated  by 

C b  -  c2/sj  ,  b  > 1/2 


VJ 


b  +  c2/sJ  ,  b  < 1/2 
,  b  -  1/2 


From:  James  M.  Dobbie/rc 
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